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(57) Abstract 

A speech signal linear prediction analysis is per- 
formed for each frame of a sp^ch signal to deteimine tl» 
coefficients of a shoit-tenn synthesis filter, and an open- 
loop analysis is performed to determine a degree of firame 
voicing. At least one closed-loop analysis is performed for 
each sub-ftame to determine an excitation sequence which, 
when applied to the short-term synthesis filter, generates a 
synthetic signal representative of the speech signal. Each 
dOGed-loop analysis uses die impulse response of a filter 
ccmsisting of the short-temi synthesis filter and a percep- 
tual weighting fito, by truncating said impulse response 
to a tnincation lengdi that is no greater tiian the number of 
samples per sub-ftame, and dependent on the energy dis- 
tributicm of said response and the degree of voicing of the 
frame. 
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On effectue, pour chaque trame de signal de parole, 
une analyse par pnSdicdon iindaire du signal de parole pour 
detemiiner les coefficients d'un filtre de synthase ft court temie, et une analyse en boucle ouvcne pour ddtcmiiner un degrtf de voisement 
de la trame. On effectue pour chaque sous-trame au moins une analyse en boucle femn^ pour ddtenniner une sequence d*excitation qui, 
soumise au filtre de synth^ k court terme, produit un signal synthdtique repr6sentatif du signal de parole. Chaque analyse en boucle 
femite utilise la r^ponse impulsionnelle d'un filtre compost du filtre de synthase k court tcrme et d*un filtre de pond^ration perceptuelle, en 
tnxiquant ladite rdponse impulsionnelle & une longueur de tjxHicature au plus dgale au nombre d*6chantillons par sous-trame et ddpendant 
de la distribution 6nerg6tique de ladite r6ponse et du degr6 de voisement de la trame. 
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PROCEDE BE CODAGE DE PAROI^E A ANAIjYSB PAR SYKTHESE 

La pr6sente invention concerne le codage de la parole 
5 utilisant 1 ' analyse par synthase. 

La demanderesse a noteonment d^crit de tels codeurs 
de parole qu'elle a d6velopp6s dans ses dexnandes de brevet 
europ6en 0 195 487, 0 347 307 et 0 469 997, 

Dans un codeur de parole k analyse par synthase, on 
10 effectue une prediction lin6aire du signal de parole pour 
obtenir les coefficients d'un filtre de synthase k court ter- 
me mod^lisant la fonction de transfert du conduit vocal. Ces 
coefficients sont transmis au d6codeur, ainsi que des param^- 
tres caract6risant une excitation k applicjuer au filtre de 
15 synthase k court terme. Dans la plupart des codeurs actuels, 
on recherche en outre les correlations ^ plus long terme du 
signal de parole pour caract6riser un filtre de synthase k 
long terme rendant compte de la hauteur tonale de la parole. 
Lorsque le signal est vois6, 1' excitation comporte en effet 

2 0 une composante predict ible pouvant etre representee par 

1 •excitation passee, retard6e de TP 6chantillons du signal 
de parole et affectee d'un gain gp. Le filtre de synthase k 
long terme, egalement reconstitue au d6codeur, a alors une 
fonction de transfert de la forme 1/B(z> avec B( z) =l-gp. z"*^^. 

25 La partie restante, non predictible, de 1' excitation est 
appeiee excitation stochastique . Dans les codeurs dits CELP 
( -Code Excited Linear Prediction" ) , 1 ' excitation stochastique 
est constituee par un vecteur recherche dans un dictionnaire 
predetermine. Dans les codeurs dits MPLPC ("Multi-Pulse 

30 Linear Prediction Coding"), 1' excitation stochastique 
comporte un certain nombre d' impulsions dont les positions 
sont recherchees par le codeur. En general, les codeurs CELP 
sont pref6res pour les bas debits de transmission, mais ils 
sont plus complexes k mettre en oeuvre que les codeurs MPLPC. 

3 5 Pour determiner le retard de prediction i long terme, 

on utilise frequeroment une analyse en boucle fermee contri- 
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buant directement k minimiser 1 • 6cart pond6r6 perceptuelle- 
ment entre le signal de parole et le signal synth^tique . 
L- inconvenient de cette analyse en boucle ferm6e est qu'elle 
est exigeante en volume de calculs, car la selection 'd' un 
retard in^slique 1' Evaluation d'un certain nombre de retards 
candidats et chaque Evaluation d'un retard necessite des cal- 
culs de produits de convolution entre 1' excitation retardEe 
et la rEponse impulsionnelle du f iltre de synthase ponderE 
perceptuellement. L • inconvenient ci-dessus existe aussi pour 
la recherche de 1" excitation stochastique. qui est egalement 
un processus en boucle fermEe oix interviennent des produits 
de convolution avec cette rEponse impulsionnelle . L' excita- 
tion varie plus rapidement que les paramEtres spectraux 
caractEristiques du f iltre de synthase h court terme. L' exci- 
tation (prEdictible et stochastique) est typiquement dEtermi- 
n6e une fois par sous-trame de 5ms. tandis que les paramfetres 
spectraux le sent une fois par trame de 20 ms. La complexity 
et la frequence de la recherche en boucle fermEe de 1 • excita- 
tion en font 1 ' Etape la plus critique quant k la rapiditE des 
calculs nEcessaires dans un codeur de parole. 

un but principal de 1' invention est de proposer un 
procEde de codage de parole de complexity rEduite en ce qui 
conceme la ou les analyses en boucle fermEe. 

L' invention propose ainsi un proc6d6 de codage h ana- 
lyse par synthase d'un signal de parole numerisE en trames 
successives subdivisEes en sous- trames comportant un nombre 
determine d« echantillons, dans lequel on effectue pour chaque 
trame une analyse par prediction lineaire du signal de parole 
pour determiner les coefficients d'un filtre de synthese k 
court terme, et une analyse en boucle ouverte pour determiner 
un degre de voisement de la trame. et on effectue pour chaque 
sous-trame au moins une analyse en boucle fermee pour deter- 
miner une sequence d' excitation qui. soumise au filtre de 
synthese h court terme, produit un signal synthetique reprE- 
sentatif du signal de parole. Chaque analyse en boucle ferm§e 
utilise la reponse impulsionnelle d'un filtre compose du fil- 
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tre de synthase k court terme et d*un filtre de pond^ration 
perceptuelle. Lors de chaque analyse en boucle fermde, on 
utilise ladite r6ponse impulsionnelle en la tronquant a une 
longueur de troncature au plus 6gale au nombre d • 6chant illons 
5 par sous-trame et dependant de la distribution energ^tique 
de ladite r6ponse et du degr6 de voisement de la trame. 

En g6n6ral , la longueur de troncature sera d'autant 
plus grande que la trame est vois6e. On peut ainsi r^duire 
sensiblement la complexity des analyses en boucle f erm^e sans 

10 perdre en quality de codage, grace ^ une adaptation aux 
caract6ristiques de voisement du signal. 

D'autres particularit6s et avantages de 1' invention 
apparaitront dans la description ci-apr^s d'exemples de 
realisation pr6f6r6s, mais non limitatifs, en r6f6rence aux 

15 dessins annexes « dans lesquels : 

- la figure 1 est un schema synoptique d'une station 
de radiocommunication incorporant un codeur de parole mettant 
en oeuvre 1* invention ; 

- la figure 2 est un schema synoptique d'une station 
20 de radiocommunication apte k recevoir un signal produit par 

celle de la figure 1 ; 

- les figures 3^6 sont des organi grammes illustrant 
un processus d' analyse LTP en boucle ouverte applique dans 
le codeur de parole de la figure 1 ; 

25 - la figure 7 est un organigramme illustrant un 

processus de determination de la r6ponse impulsionnelle du 
filtre de synthase pond6r6 appliqu6 dans le codeur de parole 
de la figure 1 ; 

- les figures 8 ^ 11 sont des organigrammes illus- 
30 trant un processus de recherche de 1' excitation stochastique 

appliqu6 da:ns le codeur de parole de la figure 1. 

Un codeur de parole mettant en oeuvre 1 ' invention est 
applicable dans divers types de syst^mes de transmission 
et/ou de stockage de parole faisant appel k une technique de 
35 con^ression niim6rique. Dans I'exemple de la figure 1, le co- 
deur de parole 16 fait partie d'une station mobile de radio- 
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communication. Le signal de parole S est un signal numdr.Que 
6chantillonn6 k une frequence typiguement 6gale h 8kHz. Le 
signal S est issu d'un convertisseur analogique-num6rxque 18 
recevant le signal de sortie amplifi6 et filtr6 d-un micro- 
5 phone 20. Le convertisseur 18 met le signal de parole S sous 
forme de trames successives elles-memes subdivis6es en-nst 
sous-trames de 1st 6chantillons . Une trame de 20 ms comporte 
typiQuement nst=4 sous-trames de lst=40 6chantillons de 16 
bits h 8kHz. En amont du codeur 16. le signal de parole S 
10 peut 6galement etre soumis k des traitements classiques de 
mise en forme tels gu'un filtrage de Hamming. Le codeur de 
parole 16 d^livre une sequence binaire de d6bit sensiblement 
plus faible que celui du signal de parole S. et adresse cette 
sequence k un codeur canal 22 dont la fonction est d' intro- 
15 duire des bits de redondance dans le signal af in de permettre 
une detection et/ou une correction d' 6ventuelles erreurs de 
transmission. Le signal de sortie du codeur canal 22 est en- 
suite moduli sur une frequence porteuse par le modulateur 24. 
et le signal modul6 est 6mis sur 1- interface axr. 
20 Le codeur de parole 16 est un codeur k analyse par 

synthase Le codeur 16 determine d'une part des param^tres 
caract6risant un f iltre de synthase h court terme mod^lisant 
le conduit vocal du locuteur, et d' autre part une sequence 
d- excitation qui, appliqu6e au filtre de synthase a court 
terme, fournit un signal synth6tique constituant une estxma-. 
tion du signal de parole S selon un crit^re de pond6ration 

perceptuelle. « -i 

A(z)=l-S3i-^ 
i=l 

Le filtre de synthase k court terme a une fonction 
de transfert de la forme 1/A(z) . avec : 

Les coefficients a^ sont determines par un module 26 
d- analyse par prediction lin^aire h court terme du signal de 
parole S. Les ai sont les coefficients de prediction lin^a^re 
35 du signal de parole S. L^ordre q de la prediction lineaire 
est typiquement de I'ordre de 10. Les methodes applicables 
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par le module 2 6 pour la prediction lin^aire A court terme 
sont bien connues dans le domaine du codage de la parole. 
Le module 26 met par exemple en oeuvre I'algorithme de 
Durbin-Levinson (voir J. Makhoul : "Linear Prediction : A 
5 tutorial review-, Proc, IEEE, Vol.63, N«>4, Avril 1975, p. 
561-580) . Les coefficients a^^ obtenus sont fournis a un 
module 28 qui les convertit en param^tres de raies spectrales 
(LSP) . La representation des coefficients de prediction a^ 
par des param^tres LSP est frequemment utilis^e dans des 

10 codeurs de parole ct analyse par synthase. Les paramditres LSP 
sont les q nombres cos(27Cf^) ranges en ordre d^croissant, les 
q frequences de raies spectrales (LSF) normalisees f^(liiiq) 
etant telles que les nombres complexes exp (27rj f ) , avec 
i = l, 3, . • . ,q-l,q+l et ^q+i = 0,5, soient les racines du 

15 polynome Q(z) defini par Q ( z ) =A ( z ) +z" '^"^"^ ^ . A ( z"^ ) et que 
les nombres complexes exp(27Cjfi), avec i = 0 , 2 , 4 , . . . , q et 
fQ=0, soient les racines du polynSme Q*(z) defini par 

Q* (2)=A(Z) -Z" l^+l> .A(2"^) . 

Les parametres LSP peuvent Stre obtenus par le module 
20 de conversion 28 par la methode classique des polynSmes de 
Cbebyshev (voir P. Kabal et R.P. Ramachandran : -The computa- 
tion of line spectral frequencies using Chebyshev polyno- 
mials-, IEEE Trans. ASSP, Vol.34, N*=> 6, 1986, pages 1419- 
1426) , Ce sont des valeurs de quantification des parametres 
25 LSP, obtenues par un module de quantification 30, qui sont 
transmises au decodeur pour que celui-ci retrouve les coeffi- 
cients a^ du filtre de synthase k court terme. Les coeffi- 
cients a^ peuvent Stre retrouves simplement, etant donne que: 

30 0(^)=(1+^"^) n (l-2cos{27cfj)z"^+z"^) 

i-1 , 3 , . • . , (j-l 



0^ (z) - d-z"^) n {1-2COS {2nf^) z'^^z'^) 

i=2,4, . • . , g 

35 
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et A{z) = [Oiz) +0*(2)]/2 

Pour 6viter des variations brusques dans la f once ion 
de transfert du filtre de synthase k court terme, les paramfe- 
tres LSP font I'objet d'une interpolation avant qu'on en de- 
duise les coefficients de prediction ai- Cette interpolation 
5 est effectu6e sur les premieres sous-trames de chaque trame 
du signal. Par exemple. si LSP^ et LSP^-i d6signent respec- 
tivement un paramfetre LSP calculi pour la trame t et pour la 
trame pr6c6dentet-l, on prend : LSP^ ( 0) =0 , 5 . LSP^.i^O . 5 . LSP^ . 
LSP^(l)=0,25.LSPt_l*0.75.LSP, et LSP^, (2 ) = . . . =LST, (nst-1) =I.SPt 

10 pour les sous-trames 0,1,2 nst-1 de la trame t. Les 

coefficients ai du filtre 1/A(z) sont alors determines, sous- 
trame par sous-trame k partir des param^tres LSP interpoles . 

Les parametres LSP non quantifies sont foumis par 
le module 28 k un module 32 de calcul des coefficients d-un 
filtre de ponderation perceptuelle 34. Le filtre de pond6- 
ration perceptuelle 34 a de pr6f6rence une fonction de 
transfert de la forme W (z) =A (z/Yi ) /A (z/72) o^I 7i et Y2 sont 
des coefficients tels que 7i>72>0 (par exemple Yi=0.9 et 
Y2=0,6). Les coefficients du filtre de ponderation percep- 
tuelle sont calcuies par le module 32 pour chaque sous-trame 
apres interpolation des paramfetres LSP regus du module 28. 

Le filtre de ponderation perceptuelle 34 re<poit le 
signal de parole S et d61ivre un signal SW pondere percep- 
tuellement qui est analyse par des modules 36. 38. 40 pour 
determiner la sequence d' excitation . La sequence d- excitation 
du filtre k court terme se compose d'une excitation 
predictible par un filtre de synthese k long terme modeiisant 
la hauteur tonale (pitch) de la parole, et d'une excitation 
stochastique non predictible, ou sequence d • innovation . 

Le module 36 effectue une prediction k long terme 
(LTP) en boucle ouverte, c'est-i-dire qu'il ne contribue pas 
directement k la minimisation de I'erreur ponderee. Dans l€ 
cas represente. le filtre de ponderation 34 intervient er 
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ainont du module d* analyse en boucle ouverte, mais il pourrait 
en etre autrement : le module 3 6 pourrait op6rer directement 
sur le signal de parole S ou encore sur le signal S 
d6barrass6 de ses correlations k court terme par un filtre 
5 de fonction de transfert A(z) . En revanche, les modules 38 
et 40 fonctionnent en boucle ferm6e, c*est-^-dire qu'ils 
contribuent directement k la minimisation de 1 ' erreur 
pond6r6e perceptuellement . 

Le filtre de synthase k long terme a une fonction de 

10 transfert de la forme 1/B(z) avec B(z) =l-gp. z"*^^ ou gp 
d6signe un gain de prediction ^ long terme et TP d6signe un 
retard de prediction ^ long terme, Le retard de prediction 
^ long terme peut typiQuement prendre N=:25 6 valeurs comprises 
entre rmin et rmax ecbantillons . Une resolution f ract ionnaire 

15 est prevue pour les plus petites valeurs de retard de facon 
k eviter les ecarts trop perceptibles en termes de frequence 
de voisement. On utilise par exemple une resolution 1/6 entre 
rmin=21 et 33 + 5/6, une resolution 1/3 entre 34 et 47^-2/3, une 
resolution 1/2 entre 48 et 88+1/2, et une resolution enti^re 

20 entre 89 et rmax=142. Chaque retard possible est ainsi 
quantifie par un index entier coir^ris entre 0 et N-l = 255. 

Le retard de prediction k long terme est determine 
en deux etapes . Dans la premiere etape, le module 3 6 d* ana- 
lyse LTP en boucle ouverte detecte les trames voisees du 

2 5 signal de parole et determine, pour chaque trame voisee, un 

degre de voisement MV et un intervalle de recherche du retard 
de prediction k long terme, Le degre de voisement MV d'une 
trame voisee peut prendre trois valeurs : 1 pour les trames 
faiblement voisees, 2 pour les trames moderement voisees, et 

3 0 3 pour les trames tr^s voisees. Dans les notations utilisees 

ci-apres, on prend un degre de voisement MV=0 pour les trames 
non voisees. L' intervalle de recherche est defini par une 
valeur centrale representee par son index de cruantif ication 
ZP et par une largeur dans le domaine des index de quantif i- 
3 5 cation, dependant du degre de voisement MV. Pour les trames 
faiblement ou moderement voisees (MV=1 ou 2 ) la largeur de 
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l-intervalle de recherche est de Nl index, c'est-a-dire que 
!■ index du retard de prediction h long terme sera recherch6 
entre ZP-16 et ZP.15 si Nl=32. Pour les trames tres vois6es 
(MV-3). la largeur de I'intervalle de recherche est de N3 
index cest-St-dire que 1' index du retard de prediction ^ 
long terme sera recherche entre ZP-8 et ZP.7 si N3=16.. 

Une fois que le degre de voisement MV d'une trame a 
6te determine par le module 36. le module 30 op^re la quanti- 
f ication des param^tres LSP qui ont auparavant ete d6termxn6s 
pour cette trame. Cette quantification est par exemple vecto- 
rielle c'est-a-dire qu'elle consiste a seiectionner , dans 
une ou plusieurs tables de quantification predeterminees . un 
jeu de parametres quantifies LSPq qui presente une distance 
minimale avec le jeu de parametres LSP fourni par le module 
28 De facon connue. les tables de quantification different 
suivant le degre de voisement MV fourni au module de quan- 
tification 30 par 1-analyseur en boucle ouverte 36. un ensem- 
ble de tables de quantification pour un degre de voisement 
MV est determine, lors d^essais prealables. de fapon ^ Stre 
statistiquement representatif de trames ayant ce degre MV. 
ces ensembles sont stockes 4 la fois dans les codeurs et dans 
les decodeurs mettant en oeuvre 1' invention. Le module 3 0 
deiivre le jeu de parametres quantifies LSPq ainsi que son 
index Q dans les tables des quantification applicables. 

Le codeur de parole 16 comprend en outre un module 
42 de calcul de la reponse in«>ulsionnelle du filtre compose 
du filtre de synthfese k court terme et du filtre de pondera- 
tion perceptuelle. Ce filtre compose a pour fonctxon d6 
transfert W(z)/A(z) . Pour le calcul de sa reponse impulsxon- 

nelle h={h(0).h(l) h(lst-l)) sur la duree d^une sous- 

trame, le module 42 prend pour le filtre de ponderatio. 
perceptuelle W(z) celui correspondant aux paramitres LSI 
interpoies mais non quantifies, c'est-^-dire celui dont le. 
coefficients ont ete calcuies par le module 32, et pour 1- 
filtre de synthese 1/A(z) celui correspondant aux paramfetre: 
LSP quantifies et interpoies. cest-^-dire celui qui ser. 
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ef f ectivement reconstitu^ par le d6codeur. 

Dans la deuxi^e 6tape de la determination du retard 
TP de prediction ^ long terme, le module 38 d' analyse LTP en 
boucle fenn6e determine le retard TP pour chaque sous-trame 
5 des trames vois6es {MV=1, 2 ou 3). Ce retard TP est 
caract6ris6 par une valeur dif f 6rent ielle DP dans le domaine 
des index de quantification, cod^e sur 5 bits si MV=1 ou 2 
(Nl = 32), et sur 4 bits si MV=3 (N3 = 16) . L ' index du retard TP 
vaut ZP+DP. De facon connue, 1* analyse LTP en boucle ferm^e 
10 consiste k determiner, dans I'intervalle de recherche des 
retards T de prediction S long terme, le retard TP qui 
maximise, pour chaque sous-trame d*une trame voisee, la 
correlation normalisee : 

15 [ J2 -«^<i) .yr*^)] 

ou x(i) designe le signal de parole ponder e SW de la sous- 

2 0 trame auquel on a soustrait la memoire du filtre de synthase 

pondere (c ' est-^-dire la reponse ct un signal nul , due ^ ses 
etats initiaux, du filtre dont la reponse impulsionnelle h 
a ete calcuiee par le module 42) , et y.p(i) designe le produit 
de convolution : 
25 i 

jf=0 

u(j-T) designant la composante predict ible de la sequence 
d* excitation retardee de T echantillons , estimee par la tech- 
30 nique bien connue du repertoire adaptatif ("adaptive 
codebook") . Pour les retards T inf6rieurs A la longueur d'une 
sous-trame, les valeurs manquantes de u(j-T) peuvent etre 
extrapoiees k partir des valeurs anterieures . Les retards 
f ractionnaires sont pris en compte en surechant illonnant le 

3 5 signal u(j-T) dans le repertoire adaptatif. Un 

surechantillonnage d*un facteur m est obtenu au moyen de 
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f litres polyphases interpolateurs . 

Le gain gp de prediction k long terme pourrait §tre 
determine par le module 38 pour chaque sous-trame. en 
appliquant la formule connue : 

J} x(i) .yj.p(i) 

_ _ ^^0 

^P' ist-1 2 



2=0 

Toutefois, dans une version pref6r6e de 1- invention, le gain 
g est calcule par le module d" analyse stochastique 40. 

L' excitation stochastique determin6e pour chaque 
sous-trame par le module 40 est de type mult i- impulsionnelle . 
une sequence d' innovation de 1st echantillons comprend np 
impulsions de positions p(n) et d'amplitude g(n). Autrement 
dit. les impulsions ont une amplitude de 1 et sont associ^es 
a des gains respectifs g(n). Etant donn4 que le retard LTP 
n-est pas determine pour les sous-trames des trames non 
vois6es. on peut prendre un nombre d' impulsions sup^rieur 
pour 1' excitation stochastique relative h ces sous-trames. 
par exemple np=5 si MV=1. 2 ou 3 et np=6 si MV=0. Les 
positions et les gains calculus par le module 40 d- analyse 
stochastique sont quantifies par un module 44. 

Un module d' ordonnancement des bits 46 recoit les 
differents param^tres qui seront utiles au decodeur, et 
constitue la sequence binaire transmise au codeur canal 22. 

Ces param^tres sont : 

- 1' index Q des paramfetres LSP quantifies pour chaque 

3 0 trame ; 

- le degre MV de voisement de chaque trame. ; 

- 1' index ZP du centre de I'intervalle de recherche 
des retards LTP pour chaque trame vois^e ; 

- 1- index differentiel DP du retard LTP pour chaque 
35 sous-trame d'une trame vois6e. et le gain associe gp ; 

- les positions p(n) et les gains g(n) des impulsions 



wo 96/21220 



PCT/FR96/000a6 



11 

de 1* excitation stochastique pour chaque sous-trame. 

Certains de ces paramfetres peuvent avoir une impor- 
tance particuli^re dans la quality de restitution de la 
parole ou une sensibility particulidre aux erreurs de 
5 transmission. On pr6voit ainsi dans le codeur un module 48 
qui recoit les dif f^rents param^tres et qui ajoute ^ certains 
d'entre eux des bits de redondance permettant de d^tecter 
et/ou de corriger d' 6ventuelles erreurs de transmission. Par 
exemple, le degr6 de voisement MV cod6 sur deux bits 6tant 

10 un paramfetre critique, on souhaite qu'il parvienne au 
d^codeur avec aussi peu d' erreurs que possible. Pour cette 
raison, des bits de redondance sont ajout6s A ce param^tre 
par le module 48. On peut par exemple ajouter un bit de 
parity aux deux bits codant MV et r6p6ter une fois les trois 

15 bits ainsi obtenus. Cet exemple de redondance permet de 
d^tecter toutes les erreurs simples ou doubles et de corriger 
toutes les erreurs simples et 75% des erreurs doubles. 

L' allocation du d6bit binaire par trame de 20 ms est 
par exemple celle indiqu6e dans le tableau I. 

20 Dans 1' exemple consid6r6 ici, le codeur canal 22 est 

celui utilise dans le systfeme paneurop6en de radiocommuni- 
cation avec les mobiles (GSM) . Ce codeur canal, d6crit en 
detail dans la Recommandation GSM 05.03, a 6te mis au point 
pour un codeur de parole ^ 13 kbit/s de type RPE-LTP qui pro- 

25 duit 6galement 260 bits par trame de 20 ms. La sensibility 
de chacun des 260 bits a 6t6 d6termin6e ^ partir de tests 
d'^coute. Les bits issus du codeur source ont 6te regroup^s 
en trois categories. La premiere de ces categories lA re- 
groupe 50 bits qui sont cod6s convolutionnellement sur la 

3 0 base d'un polyn6me g6n6rateur donnant une redondance d'un 
demi avec une longueur de contrainte 6gale k 5. Trois bits 
de parite sont calculus et ajout^s aux 50 bits de la cat6- 
gorie lA avant le codage convolutionnel . La seconde cat^gorie 
(IB) conpte 132 bits qui sont prot6g6s ^ un taux d'un demi 

3 5 par le m§me polyndme que la cat6gorie precedent e. La 
troisi^me cat6gorie (ID contient 7 8 bits non proteges. Apr^s 
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application du code convolutionnel , les bits (456 par trame) 
sont soumis ^ un entrelacement . Le module d- ordonnancement 
46 du nouveau codeur source mettant en oeuvre 1' invention 
distribue les bits dans les trois categories en fonctio:.. de 
1 ■ importance subjective de ces bits. 
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ZP 
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DP 
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24 1 


1 positions 


80 




72 


1 72 1 


1 impulsions 










1 gains 


140 




100 


100 1 


1 impulsions 










1 Total 


260 


260 


I 260 1 



une station mobile de radiocommunication apte a rece- 
voir le signal de parole traits par le codeur source 16 est 
representee schema tiquement sur la figure 2. Le signal radio 
recu est d-abord traite par un d^modulateur 50 puis par un 
decodeur canal 52 qui effectuent les operations duales de 
celles du modulateur 24 et du codeur canal 22. Le decodeur 
canal 52 fournit au decodeur de parole 54 une sequence b.- 
naire qui. en 1' absence d'erreurs de transmission ou lorsque 
les eventuelles erreurs ont ^t^ corrig6es par le decodeur 
canal 52, correspond a la sequence binaire qu'a d41ivree le 
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module d' ordonnancement 46 au niveau du codeur 16. Le d6co- 
deur 54 comprend un module 56 qui reqoit cette sequence bi- 
naire et qui identifie las param^tres relatifs aux dif f^ren- 
tes trames et sous-trames. Le module 56 effectue en outre 
5 quelques contrSles sur les paramfetres recus . En particulier, 
le module 56 examine les bits de redondance introduits par 
le module 48 du codeur, pour d6tecter et/ou corriger les 
erreurs affectant les param^tres associ6s ^ ces bits de 
redondance. 

10 Pour chaque trame de parole ^ synth6tiser, un module 

58 du d^codeur recoit le degr6 de voisement MV et 1 ' index de 
Q de quantification des param^tres LSP. Le module 58 retrouve 
les paramfetres LSP quantifies dans les tables correspondant 
A la valeur de MV, et , aprfes interpolation, les convertit en 

15 coefficients a^ pour le filtre de synthase ^ court terme 60. 

Pour chaque sous- trame de parole k synth6tiser, un g6n6rateur 
d' impulsions 62 regoit les positions p{n) des np impulsions 
de !• excitation stocbastique. Le g6n6rateur 62 d61ivre des 
impulsions d' amplitude unitaire qui sont cbacune multipli6es 

20 en 64 par le gain associ6 g(n) . La sortie de 1 ' amplif icateur 
64 est adressee au filtre de synthese ^ long terme 66. Ce 
filtre 66 a une structure h repertoire adaptatif, Les 
echantillons u de sortie du filtre 66 sont m4moris§s dans le 
repertoire adaptatif 68 de facon h §tre disponibles pour les 

25 sous-trcones ult6rieures. Le retard TP relatif i une sous- 
trame, calculi i partir des index de quantification ZP et DP, 
est fourni au repertoire adaptatif 68 pour produire le signal 
u convenablement retarde. L ' amplif icateur 70 multiplie le 
signal ainsi retarde par le gain gp de prediction ^ long 

30 terme. Le filtre k long terme 66 comprend enfin un 
additionneur 72 qui ajoute les sorties des amplif icateurs 64 
et 70 pour fournir la sequence d'excitation u. Lorsque 1 ' ana- 
lyse LTP n'a pas ete effectuee au codeur, par exemple si 
MV=0, un gain de prediction gp nul est impose a 1 'amplif i- 

35 cateur 70 pour les sous-trames correspondant es . La sequence 
d'excitation est adressee au filtre de synthase ^ court terme 
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60, et le signal resultant peut encore, de facon connue, Stre 
somnis k un posc-filtre 74 dont les coefficients dependent 
des paramfetres de synthase recus, pour former le signal de 
parole synth^tique S • . Le signal de sortie S' du d6codedr 54 
5 est ensuite convert! en analogique par le convert isseur 7 6 
avant d'Stre amplifi^ pour commander un haut-parleur 78. 

on va maintenant d6crire, en reference au figures 3 
h 6, le processus d' analyse LTP en boucle ouverte mis en 
oeuvre par le module 3 6 du codeur suivant un premier aspect 

10 de 1' invention. 

Dans une premiere etape 90, le module 36 calcule et 

memorise, pour chaque sous-trame st=0,l nst-1 de la trame 

courante, les autocorrelations Cg^Oc) et les Energies retar- 
d6es Gst(k) du signal de parole pondere SW pour les retards 

15 entiers k compris entre rmin et rmax : 

(st+1) .lst-1 
C Ak)= S .SW{i-k) 
^ i'St.lst 

(sc+l).ist-l 

20 G^Ak)^ E [SW(i-k)t 

i=st.lat 

Les Energies par sous-trame RO^^ sont 6galement 

calcul6es : . , ^ , 

(St+1) .IsC-1 

i=at.lst 

A l'6tape 90, le module 36 determine en outre, pour 
chaque sous-trame st, le retard entier Kg^ Qui maximise 
1- estimation en boucle ouverte P^^i)^) du gain de prediction 
^ long terme sur la sous-trame st, en excluant les retards 
k pour lesquels 1 • autocorrelation C^^W est negative ou plus 
petite qu'une petite fraction e de 1 • 6nergie RO^. de la sous- 



30 



trame. L-estimation Pst«^' exprimee en ddcibels s ' 6crit : 
Maximiser P^^i)^) revient done h maximiser I'expressior 
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XQ^(k)=C^^^ (k) /GQf^(k) comme indiqu6 sur la figure 6. Le 
retard entier Kg^. est le retard de base en resolution enti^re 
pour la sous-trame st. L'6tape 90 est suivie par une 
comparaison 92 entre une premiere estimation en boucle 
5 ouverte du gain de prediction global sur la trame courante 
et un seuil predetermine SO typiquement compris entre 1 et 
2 decibels (par exemple SO =1,5 dB) . La premiere estimation 
du gain de prediction global est egale k : 

10 r ^^^'^ 1 

8t=0 

ou RO est I'energie totale de la trame 

(RO =: ROo+ ROi+.,.+ ROnst-l^' ^st^^st^^^st^^'^st^/^st^^'st:^ 

15 designe le maximum determine k I'etape 90 relativement ci la 
sous-trame St. Comme 1 ' indique la figure 6, la comparaison 
92 peut etre ef fectuee sans avoir k calculer le logarithme. 
Si la comparaison 92 montre une premiere estimation 
du gain de prediction inferieure au seuil SO, on considere 

20 que le signal de parole contient trop peu de correlations k 
long terme pour etre voise, et le degre de voisement MV de 
la trame courante est pris egal k 0 k 1 ' etape 94, ce qui 
termine dans ce cas les operations effectuees par le module 
36 sur cette trame • Si au contraire le seuil SO est depass6 

25 k 1 ' etape 92, la trcime courante est detectee comme voisee et 
le degre MV sera egal ^1* 2 ou 3 . Le module 36 calcule 
alors, pour chaque sous-trame st, une liste Ig^. contenant des 
retards candidats pour constituer le centre ZP de 1 ' inter- 
valle de recherche des retards de prediction k long terme. 

3 0 Les operations effectu6es par le module 3 6 pour 

chaque sous-trame st (st initialise ct 0 ^ 1 ' etape 96) d'une 
trame voisee comroencent par la determination 98 d'un seuil 
de selection SEg^ en decibels egal k une fraction determinee 
P de 1' estimation ^st^^st^ gain de prediction en decibels 

35 sur la sous-trame, maximis6e k l*6tape 90 (p=0,75 typique- 
ment). Pour chaque sous-trame st d'une trame voisee, le 
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module 36 determine le retard de base rbf en resolution 
enti^re pour la suite du traitement. Ce retard de base 
pourrait Stre pris 6gal ^ I'entier K^^ obtenu ^ 1 • 4tape 90. 
Le fait de rechercher le retard de base en resolution 
fractionnaire autour de K^^ permet toutefois de gagner en 
precision. L-etape 100 consiste ainsi. ^ rechercher. autrour 
du retard entier K^^ obtenu ^ I'^tape 90. le retard 
fractionnaire qui maximise 1 • expression C^t /Ggf ^ette 
recherche peut gtre effectu6e h la resolution maximale des 
retards f ractionnaires (1/6 dans 1 ■ exemple decrit ici) m§me 
si le retard entier K^^ n-est pas dans le domaine ou cette 
resolution maximale s' applique. On determine par exemple le 
nombre A^, qui maximise C^^^ (K,,.b/6) /G,,(K,,.S/6) pour - 
6<5<.6. puis le retard de base rbf en resolution maximale est 
pris egal k K^^* ^st^'^" valeurs f ractionnaires T du 

retard, les autocorrelations C^^iT) et les energies retardees 
G (T) sont obtenues par interpolation ^ partir des valeurs 
memorisees ^ 1 " etape 90 pour les retards entiers. Bien 
entendu, le retard de base relatif ^ une sous-trame pourrait 
egalement gtre determine en resolution fractionnaire des 
l-etape 90 et pris en compte dans la premiere estimation du 
gain de prediction global sur la trame. 

una fois que le retard de base rbf a ete determine 
pour une sous-trame. on precede A un examen 101 des sous- 
multiples de ce retard af in de retenir ceux pour lesquels le 
gain de prediction est relativement important (figure 4). 
puis des multiples du plus petit sous-multiple retenu (figure 
5) A 1-etape 102. 1 • adresse j dans la liste 1^^ et 1 • index 
mdu sous-multiple sont initialises ^ 0 et 1. respect ivement . 
une comparaison 104 est effectuee entre le sous-multiple 
rbf/m et le retard minimal rmin. Le sous-multiple rbf/m est 
^ examiner s'il est superieur ^ rmin. On prend alors pour 
I'entier i la valeur de 1 • index du retard quantifie r^ le 
plus proche de rbf/m (etape 106), puis on compare, en 108, 
la valeur estimee du gain de prediction P^^ir^) associee au 
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retard quantifi6 r^^ pour la sous-trame consid6r6e au seuil 
de selection SEg^ calcul6 ^ l'6tape 98 : 

Pst(ri) = 20.1ogj^Q[RO^^/[RO^^' C^^^ (r ^) /G^^ (r ] ] 
avec, pour les retards f ractionnaires une interpolation des 
5 valeurs Cg^^ et Gg^. calcul6es ^ 1 * 6tape 90 pour les retards 
entiers. Si Pg^. (r^^) <SE3^ , le retard r^ n'est pas pris en 
consideration, et on passe directement ^ 1 • 6tape 110 
d' incr6inentation de 1 ' index m avant d'effectuer de nouveau 
la comparaison 104 pour le sous-multiple suivant . Si le test 

10 108 montre que Pg^(r^) i SEg^., le retard r^ est retenu et on 
execute l'6tape 112 avant d ' incr^menter 1 • index m h 1 ' 6tape 
110. A 1 • 6tape 112, on memorise 1 ' index i A I'adresse j dans 
la liste Igf on donne la valeur m k I'entier mO destine ^ 
etre 6gal ci 1 ' index du plus petit sous-multiple retenu, puis 

15 on incr^mente d*une unit^ I'adresse j. 

L'examen des sous-multiples du retard de base est 
termini lorsque la comparaison 104 montre rbf/m < rmin. On 
examine alors les retards multiples du plus petit rbf/mO des 
sous-multiples pr^c^demment retenus suivant le processus 

20 illustr6 sur la figure 5. Cet examen commence par une initia- 

lisation 114 de 1* index n du multiple : n=2 . Une comparaison 
116 est effectu6e entre le multiple n.rbf/mO et le retard 
maximal rmax. Si n.rbf/mO > rmax, on effectue le test 118 
pour determiner si 1 ' index mO du plus petit sous-multiple est 

25 un multiple entier de n. Dans 1 ' affirmative, le retard 
n.rbf/mO a d6j^ 6te examine lors de 1 • examen des sous- 
multiples de rbf, et on passe directement k 1 • etape 12 0 
d ' incrementation de 1 ' index n avant d'effectuer de nouveau 
la comparaison 116 pour le multiple suivant. Si le test 118 

30 montre que mO n'est pas un multiple entier de n, le multiple 
n.rbf/mO est i examiner. On prend alors pour 1' entier i la 
valeur de 1 ' index du retard quantifie r^^ le plus proche de 
n.rbf/mO {6tape 122), puis on compare, en 124, la valeur 
estimee du gain de prediction Pg^(rj^) au seuil de selection 

35 ^^sf ^st ^^i^'^^^st ' ^® retard r^^ n'est pas pris en consi- 
deration, et on passe directement ^ l'6tape 120 d'incremen- 
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tation de 1 ' index n. Si le test 124 montre que 
p (r - ) i SEgt' retard r^ est retenu et on execute I'^tape 
126 atant d' increment er 1 • index n k l'6tape 120. A I'^tape 
126, on memorise 1 ' index i h I'adresse j dans la lisce 1^^. 
puis on incr^mente d'une unit6 1 ■ adresse j. 

L'examen des multiples du plus petit sous-multiple 
est termini lorsque la comparaison 116 montre que 
n.rbf/mO > ntiax. A ce moment, la liste Ig^ contient j index 
de retards candidats. Si on souhaite limiter a jmax la lon- 
gueur maximale de la liste Igt PO^^ Stapes suivantes, on 
peut prendre la longueur jg^ ^ette liste 6gale a 
min(j.jmax) (6tape 128) puis, a 1 • 6tape 130, ordonner la 
liste «ians I'ordre des gains C^t^ (rj^C (j ) ^ '^^St (^lst(j)^ 
d^croissants pour Osj<jst facon a ne conserver que les jg^ 
retards procurant les plus grandes valeurs de gain. La valeur 
de jmax est choisie en fonction du compromis vis6 entre 1-ef- 
f icacit6 de la recherche des retards LTP et la complexity de 
cette recherche. Des valeurs typiques de jmax vont de 3 ^ 5 . 

Une fois que les sous-multiples et les multiples ont 
6t6 examines et que la liste Ig^ ^ ainsi 6t6 obtenue (figure 
3) , le module d' analyse 3 6 calcule une quantity Ymax determi- 
nant une seconde estimation en, boucle ouverte du gain de 
prediction k long terme sur 1- ensemble de la trame, ainsi que 
des index ZP, ZPO et ZPl dans une phase 132 dont le d6roule- 
ment est d6taill6 sur la figure 6. Cette phase 132 consiste 
a tester des intervalles de recherche de longueur Nl pour 
determiner celui qui maximise une deuxi^e estimation du gain 
de prediction global sur la trame. Les intervalles testes 
sont ceux dont les centres sont les retards candidats conte- 
nus dans la liste Igt calculee lors de la phase 101. La phase 
132 commence par une etape 136 ou 1 • adresse j dans la liste 
I ^ est initialisee ^ 0. A 1 ' etape 138, on v6rif ie si 1 • index 
(j) a dej^ 6te rencontre en testant un intervalle prece- 
dent centre sur Ist'<3') avec sf<st et Osj'<jsf' ^^^^^ 
d-eviter de tester deux fois le mSme intervalle. Si le test 
138 reveie que Igt^^^ figurait d^^k dans une liste 1^^. avec 
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st'<st, on incr^mente directement I'adresse j ^ 1 ' 6tape 140, 
puis on la compare ^ la longueur jg^ de la liste Si la 

comparaison 142 montre que j<jg^., on revient k 1 ' 6tape 138 
pour la nouvelle valeur de I'adresse j. Lorsque la comparai- 

5 son 142 montre que j^jg^-, tous les intervalles relatifs k la 
liste Ig^. ont 6t6 testes, et la phase 132 est termin6e. Lors- 
quB le test 138 est n6gatif, on teste 1 ' inteirvalle centr6 sur 
Ig^{j) en commencant par l'6tape 148 ou on determine, pour 
chaque sous-trame st • , 1 ' index igt ' retard optimal qui 

0 maximise sur cet intervalle 1' estimation en boucle ouverte 
Pg^(rj^) du ,gain de prediction ^ long terme, c'est-^- 

dire qui maximise la quantity ^st*^^^~^st*^^^i^^^st'^^i^ 
ou r^^ d^signe le retard quantifi6 d* index i pour 
Ig^.{j)-Nl/2 ii<Ig^ ( j ) +N1/2 et Oii<N. Lors de la maximisation 
5 148 relative ^ une sous-trame st ' , on 6carte a priori les 
index i pour lesquels 1 • autocorrelation Cg^.(r^> est nega- 
tive, pour 6viter de d6grader le codage. S'il se trouve que 
toutes les valeurs de i comprises dans 1' intervalle teste 
[I(j)-Nl/2, I(j)+Nl/2[ donnent lieu ^ des autocorrelations 
20 Cg^..(r^) negatives, on seiectionne 1 ' index ig^-. pour lequel 
cette autocorrelation est la plus petite en valeur absolue. 
Ensuite, en 150, la quantite Y determinant la deuxifeme esti- 
mation du gain de prediction global pour 1' intervalle centre 
sur lgt;tj) ®st calcuiee selon : 

jist-1 

puis compar6e ^ Ymax, oCi Ymax represent e la valeur ^ 
30 maximiser. Cette valeur Ymax est par exemple initialisee ^ 
0 en m&me temps cjue 1 • index st ci I'etape 96. Si YiYmax, on 
passe directement cl I'etape 140 d' incrementation de 1 • index 
j . Si la con^araison 150 montre que Y > Ymax, on execute 
I'etape 152 avant d' incrementer I'adresse j ^ I'etape 140, 
35 A cette etape 152, 1 ' index ZP est pris egal i Ig^(j) et les 
index 2P0 et ZPl sont respect ivement pris egaux au plus petit 



PCT/FR96/00006 

y/O 96/21220 



20 



10 



15 



et au plus grand des index d^tenain^s I'^tape 148. 

A la fin de la phase 132 relative k une sous-trame 
St 1' index St est incr6inent6 d'une unit6 (6tape 154) puis 
compar6. h 1 • 6tape 156, au nombre nst de sous- trames - par 
trame. Si st<nst, on revient k I'^tape 98 pour effectuer les 
operations relatives k la sous-trame suivante. Lorsque la 
coiaparaison 156 montre que st=nst, i; index ZP d^signe le 
centre de I'intervalle de recherche qui sera fourni au module 
38 d' analyse LTP en boucle ferm6e, et ZPO et ZPl sont des 
index dont 1 ' 6cart est repr6sentatif de la dispersion des 
retards optimaux par sous-trame dans 1 • intervalle centr6 sur 
ZP . 

A I'^tape 158, le module 36 determine le degr6 de 
voisement MV, sur la base de la seconde estimation en boucle 
ouverte du gain exprim^e en decibels : Gp=20. log^Q (RO/^O- 
Ymax) . on fait appel k deux autres seuils SI et S2 . Si GpsSl. 
le degr^ de voisement MV est pris egal ^ 1 pour la trame cou- 
rante. Le seuil SI est typiquement compris entre 3 et 5 dB; 
par exemple Sl=4 dB. Si Sl<Gp<S2, le degr6 de voisement MV 
20 est pris 6gal ^ 2 pour la trame courante. Le seuil S2 est ty- 
piquement compris entre 5 et 8 dB ; par exemple S2=7 dB. Si 
Gp>S2. on examine la dispersion des retards optimaux pour les 
diff6rentes sous-trames de la trame courante. Si ZPl-ZP<N3/2 
et ZP-ZP0SN3/2. un intervalle de longueur N3 centr6 sur ZP 
suffit a prendre en compte tous les retards optimaux et le 
degr6 de voisement est pris 6gal k 3 (si Gp>S2) . Sinon. si 
ZPl-ZPiN3/2 ou ZP-ZPO>N3/2. le degr6 de voisement est pris 

6gal a 2 (si Gp>S2) . 

L- index ZP du centre de 1' intervalle de recherche du 
retard de prediction pour une trame vois6e peut etre compris 
entre 0 et N-l=255, et 1 ' index diff^rentiel DP determine pour 
le module 38 peut aller de -16 k .15 si MV=1 ou 2 . et de -8 
k .7 Si MV=3 (cas Nl=32. N3 = 16) . L' index ZP-.DP du retard TP 
finalement d6termin6 peut done dans certains cas itre plus 
35 petit que 0 ou plus grand que 255. Ceci permet k 1' analyse 
LTP en boucle f erm6e de porter 6galement sur quelques retards 
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TP plus pet its que rmin ou plus grands que rmax. On ain61iore 
ainsi la quality subjective de la restitution des voix dites 
pathologiques at des signaux non vocaux (frequences vocales 
DTMF ou frequences de signalisation utilis6es par le reseau 
5 teiephonique commut6) . Une autre possibility est de prendre 
pour I'intervalle de recherche les 32 premiers ou derniers 
index de quantification des retards si ZP<16 ou ZP>240 avec 
MV-1 ou 2 / et les 16 premiers ou derniers index si ZP<8 ou 
ZP>248 avec MV=3. 
10 Le fait de r6duire I'intervalle de recherche des 

retards pour les trames tr^s vois^es (typiquement 16 valeurs 
pour MV=3 au lieu de 32 pour MV=1 ou 2) permet de diminuer 
la complexity de 1 ' analyse LTP en boucle ferm6e effectu^e par 
le module 38 en r6duisant le nombre de convolutions y»p(i) a 
15 calculer suivant la formule (1). Un autre avantage est qu ' un 
bit de codage de 1 ' index diff6rentiel DP est economist- Le 
debit de sortie etant constant, ce bit peut §tre realloue au 
codage d'autres pararndtres. On peut en particulier allouer 
ce bit suppiementaire k la quantification du gain de pr^dic- 
20 tion k long terme gp calcuie par le module 40. En effet, une 
meilleure precision sur le gain gp grSce ^ un bit de quanti- 
fication suppiementaire est appreciable car ce parametre est 
perceptuellement important pour les sous-trames tr^s voisees 
(MV=3) . Une autre possibilite est de pr6voir un bit de parite 
25 pour le retard TP et/ou le gain gp, permet tant de detecter 
d • eventuelles erreurs affectant ces parametres. 

11 est possible d'apporter quelques modifications 
au processus d* analyse LTP en boucle ouverte decrit ci-dessus 
en reference aux figures 3^6, 
30 Suivant une premiere variante de ce processus, les 

premieres optimisations effectuees ^ I'etape 90 relativement 
aux diff6rentes sous-trames sont remplacees par une seule 
optimisation portant sur 1* ensemble de la trame. Outre les 
parametres Cg^^Oc) et Gg^<k) calcuies pour chaque sous-trame 
35 St, on calcule egalement les autocorrelations C{k) et les 
energies retardees G{k) pour 1* ensemble de la trame : 
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JlSt-1 

c(Jc)= 5: c^^{k) 
st=0 
nst-l 

5 G{k)= S G^^(JC) 

st=0 

On determine alors le retard de base en resolution 
entifere K qui maximise X (k) =C^ (k) /G (k) pour rmin i k S rmax. 
10 La premiere estimation du gain comparee a SO ^ 1 • 6tape 92 est 
alors p{K)=20.1ogio^^°^ "^^^"^^^^ ^ ^ • On determine ensuite, au- 
tour de K, un seul retard de base en resolution f ractionnaire 
rbf et I'examen 101 des sous-multiples et des multiples est 
effectue une seule fois et produit une seule liste I au lieu 
15 de nst listes Igf ^a phase 132 est ensuite effectu6e une 
seule fois pour cette liste I, en ne distinguant les sous- 
trames qu'aux 6tapes 148. 150 et 152. Cette variante de 
realisation a pour avantage de r^duire la complexite de 
1" analyse en boucle ouverte. 
20 Suivant une seconde variante du processus d* analyse 

LTP en boucle ouverte. le domaine I rmin, rmax] des retards 
possibles est subdivis6 en nz sous-intervalles ayant par 
exemple la meme longueur (n2=3 typiquement) . et les premieres 
optimisations effectu^es h 1 • etape 90 relativement aux diffe- 
25 rentes sous-trames sont remplac^es par nz optimisations dans 
les differents sous-intervalles portant chacune sur !• ensem- 
ble de la trame. On obtient ainsi nz retards de base 
V • K_ ' en resolution entiere. La decision voise/non 

voise (etape 92) est prise sur la base de celui des retards 
3 0 de base K^' qui procure la plus grande valeur pour la 
premiere estimation en boucle ouverte du gain de prediction 
k long terme. Ensuite. si la trame est voisee, on determine 
les retards de base en resolution f ractionnaire par le mSme 
processus qu'^ 1 • etape 100, mais en autorisant seulement les 
35 valeurs de retard quantifiees. L'examen 101 des sous- 
multiples et des multiples n'est pas effectue. Pour la phase 



wo 96/21220 



PCT/FR96/00006 



23 

132 de calcul de la seconde estimation du gain de prediction, 
on prend comme retards candidats les nz retards de base pr6- 
c6deininent determines. Cette seconde variante permet de se 
dispenser de 1 ' examen syst6matique des sous-multiples et des 
5 multiples qui sont en general pris en consideration grace ^ 
la subdivision du domaine des retards possibles. 

Suivant une troisi^e variante du processus d' analyse 
LTP en boucle ouverte, la phase 132 est modifiee en 
ce que, aux etapes d* optimisation 148, on determine 

10 d'une part 1 • index ig^. , qui maximise (r^^) /G^^ . (r^^) pour 

Ist<^^'^^^^^^<^st(^'^^^^^2 0^i<N, et d» autre part, au 

cours de la meme boucle de maximisation, 1 • index k^^, qui 
maximise cette meme quantite sur un intervalle reduit 
Igt <3>"N3/2ii<I^^rj;+W3/2 et Oii<N, L^etape 152 est 

15 egalement modifi6e : on ne memorise plus les index ZPO et 
ZPl, mais une quantite Ymax ' d6finie de la mSme mani^re que 
Ymax mais en reference ^ 1 ' intervalle de longueur reduite : 

nst-1 

20 ^niax'= E ^sc'^*st'^ 

Dans cette troisieme variante, la determination 158 
du mode de voisement conduit h select ionner plus souvent le 
degre de voisement MV=3. On prend egalement en compte, en 

2 5 plus du gain Gp precedemment decrit, une troisieme estimation 
en boucle ouverte du gain LTP, correspondant k Ymax' : 
Gp'=20.1og-|^Q(R0/(R0-Ymax' ) 1 . Le degre de voisement est MV=1 
si GpstSl, MV=3 si Gp'>S2 et MV=2 si aucune de ces deux 
conditions n'est verifiee. En augmentant ainsi la proportion 

30 de trames de degre MV=3, on reduit la complexite moyenne de 
1* analyse en boucle fermee et on ameiiore la robustesse aux 
erreurs de transmission, 

Une quatri^e variante du processus d* analyse LTP en 
boucle ouverte concerne surtout les trames faiblement voisees 

35 (MV=l) , Ces trames correspondent souvent 3i un debut ou ^ une 
fin d'une zone de voisement. Frequemment , ces trames peuvent 
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coinporter de une h trois sous-trames pour lesquelles le coef- 
ficient de gain du filtre de synthase ^ long terme est nul 
voire n6gatif . II est propose de ne pas effectuer 1' analyse 
LTP en boucle fermee pour les sous-trames en question, afin 
5 de r6duire la complexity moyenne du codage. Ceci peut atre 
r6alis6 en m^morisant k 1 • 6tape 152 de la figure 6 nst poin- 
teurs indignant pour chaque sous-trame st • si 1 • autocorrela- 
tion Cgt' correspondant au retard d' index L^^. est negative 
ou encore tr^s petite. Une fois que tous les interval les r6- 
10 f6renc6s dans les listes 1^^. les sous-trames pour lesquelles 
le gain de prediction est n6gatif ou n^gligeable peuvent Stre 
identifiees en consulcant les nst pointeurs. Le cas echeant 
le module 3 8 est d^sactive pour les sous-trames correspon- 
dantes. Ceci n'affecte pas la quality de 1 • analyse LTP puis- 
15 que le gain de prediction correspondant ^ ces sous-trames 
sera de toutes f aeons quasiment nul . 

Un autre aspect de 1' invention concerne le module 42 
de calcul de la rdponse impulsionnelle du filtre de synthese 
20 pondere. Le module 38 d' analyse LTP en boucle fermee a besoin 
de cette r^ponse impulsionnelle h sur la dur^e d'une sous- 
trame pour calculer les convolutions y^Ci) selon la f ormule 
(1). Le module 40 d'analyse stochastique en a egalement 
besoin pour calculer des convolutions comme on le verra plus 
25 loin. Le fait d' avoir k calculer des convolutions avec une 
reponse h s'^tendant sur la duree d'une sous-trame (lst=40 
typiquement) implique une relative complexity du codage. 
qu'il serait souhaitable de r^duire notamment pour augmenter 
1-autonomie de la station mobile. Dans certains cas il a ete 
30 propose de tronquer la reponse impulsionnelle h une longueur 
inferieure h la longueur d'une sous-trame (par exemple ^ 20 
echantillons) , mais ceci peut degrader la qualit6 du codage. 
on propose selon 1' invention de tronquer la reponse impul- 
sionnelle h en tenant compte d'une part de la distributior 
35 energetique de cette reponse et d' autre part du degre de voi- 
sement MV de la trame consideree. determine par le module 3€ 
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d' analyse LTP en boucle ouverte. 

Les operations ef features par le module 42 sont par 
exemple conformes k 1 ' organigramme de la figure 7. La r^ponse 
impulsionnelle est d'abord calcul6e ^ 1 ' 6tape 160 sur une 
5 longueur pst sup^rieure ^ la longueur d'une sous-trame et 
suf f isainment grande pour qu'on soit assur6 de prendre en com- 
pte toute l*6nergie de la r^ponse impulsionnelle (par exemple 
pst = 60 pour nst=4 et lst = 40 si la prediction lin^aire ^ court 
terme est d'ordre q=10) . A I'^tape 160, on calcule ^galement 
10 les energies tronqu^es de la r^ponse impulsionnelle : 

i 
ic=0 

Les composantes h(i> de la r^ponse impulsionnelle et 
15 les energies tronquees Eh(i) peuvent &tre obtenues en 
filtrant une impulsion unitaire au moyen d'un filtre de 
fonction de transfert W(z)/A(2) d*etats initiaux nuls, ou 
encore par recurrence : 



20 
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f(i)=6(i)+£ a^[Y*.f(i-ic)- yf-SCi-Jt)] (2) 

Q 

h{i) = fU) + 2 ^k'^^^^'^^ (3) 
Eb(i) « sh(l'l) ^ lh(i)]^ 



pour 0<i<pst, avec f{i)=h(i)=0 pour i<0, 5 { 0 ) =f ( 0 ) =h ( 0 ) 
=Eh{0)=l, et 5(i)=0 pour i#0. Dans 1 'expression (2), les 
coefficients aj^ sont ceux intervenant dans le filtre de 
ponderation perceptuelle, c'est-i-dire les coefficients de 

30 prediction lineaire interpoies mais non quantifies, tandis 
que dans 1' expression (3), les coefficients aj^ sont ceux 
appliques au filtre de synthese, c'est-a-dire les 
coefficients de prediction lineaire quantifies et interpoies. 

Ensuite le module 42 determine la plus petite lon- 

35 gueur La telle que I'energie Eh(La-l) de la reponse impul- 
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sionnelie tronqu6e ^ La 6chantillons soit au moins 6gale ^ 
une proportion a de son 6nergie totals Eh{pst-1) estim^e sur 
pst 6chantillons. Une valeur typique de a est 98%. Le noxnbre 
La est initialise k pst h l-6tape 162 et decrement^ d'une 
5 unite en 166 tant que Eh (La-2 ) >a. Eh (pst-1) (test 164).. La 
longueur La cherch^e est obtenue lorsque le test 164 montre 
que Eh(La-2)<a.Eh(pst-i) . 

Pour tenir compte du degr^s de voisement MV, un terme 
correcteur A{MV) est ajout6 h la valeur de La qui a 6t6 
10 obtenue (^tape 168) . Ce terme correcteur est de preference 
une fonction croissante du degre de voisement. On peut par 
exemple prendre A(0) =-5. A ( 1 ) =0 . A { 2 ) = + 5 et A{3)=-.7. De cette 
fagon. la reponse impul sionnelie h sera determinee de facon 
d'autant plus precise que le voisement de la parole est 
15 important. La longueur de troncature Lh de la reponse impul- 
sionnelle est prise egale 4 La si Lasnst et k nst sinon. Les 
echantillons restants de la reponse impulsionnelle (h(i)=0 
avec iiLh) peuvent dtre annuies . 

Avec la troncature de la reponse impulsionnelle, le 
2 0 calcul (1) des convolutions yr^li) par le module 38 d' analyse 
LTP en boucle fermee est modi fie de la facon suivante : 

:f=inax(0,i-Lh+l) 

L'obtention de ces convolutions, qui represente uh€ 
part important e des calculs effectues, necessite done sensi- 
blement moins de mult iplications , d'additions et d'adressages 
dans le repertoire adaptatif lorsque la reponse impulsion- 
nelle est tronquee. La troncature dynamique de la reponse im- 
pulsionnelle faisant intervenir le degre de voisement MV per- 
met d'obtenir une telle reduction de complexite sans affecte: 
la qualite du codage. Les m§mes considerations s'appliquen 
pour les calculs de convolutions effectues par le module 4 
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d- analyse stochastique. Ces avantages sont particuli^rement 
appr6ciables lorsque le f litre de pond^ration perceptuelle 
a une fonction de transfert de la forme Vt (z) =K(,z/y^) /k(z/-^2^ 
avec 0<Y2<7i<l qui donne lieu des r^ponses impulsionnelles 
5 g6n6ralement plus longues que celles de la forme 
W(z)=A(z) /A(2/Y> plus commun^ment employees dans les codeurs 
h analyse par synthase. 

Un troisi^me aspect de 1* invention concerne le 

10 module 40 d' analyse stochastique servant k modeliser la 
partie non pr^dictible de 1 ' excitation . 

L" excitation stochastique consid§r6e ici est de type 
multi-impulsionnelle. L* excitation stochastique relative k 
une sous-trame est representee par np impulsions de positions 

15 p(n) et d'amplitudes, ou gains, g(n) (isnsnp) . Le gain gp de 
prediction k long terme peut 6galement §tre calcul6 au cpurs 
du mSme processus. De facon g^n^rale, on peut consid^rer que 
la sequence d' excitation relative k une sous-trame comporte 
nc contributions associ^es respect ivement k nc gains. Les 

20 contributions sont des vecteurs 1st echantillons qui, pon- 
d6r6s par les gains associes et somm^s correspondent k la 
sequence d' excitation du filtre de synthase k court terme. 
Une des contributions peut Stre predictible, ou plusieurs 
dans le cas d'un filtre de synthase k long terme k plusieurs 

25 prises ( -multi- tap pitch synthesis filter-) . Les autres con- 
tributions sont dans le cas present np vecteurs ne comportant 
que des 0 sauf une impulsion d' amplitude 1 . On a done nc=np 
si MV=0, et nc=np+l si MV=1, 2 ou 3 . 

L' analyse multi-impulsionnelle incluant le calcul du 

3 0 gain gp=g(0) consiste, de facon connue, k trouver pour chaque 
sous-trame des positions p(n) (isninp) et des gains g(n) 
(Osninp) qui minimisent 1 • erreur quadratique ponderee 
perceptuellement E entre le signal de parole et le signal 
synthetase, donnee par : 

35 
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nc-1 .2 

les gains 6tant solution du systfeme lin^aire g.B=b. 
Dans les notations ci-dessus : 

- X d^signe un vecteur-cible initial compost des 1st 
6chantillons du signal de parole pond6r6 SW sans m^moire : 

X= (x(0) ,x(l) x(lst-l) ) , les x(i) ayant 6t6 calculus cornme 

indique pr#c6deinment lors de 1' analyse LTP en boucle ferm^e ; 

- g d§signe le vecteur ligne compost des np+1 gains : 

g=(g(0)=gp, gd) g(np)) ; 

- les vecteurs-ligne ^^in) (OSn<nc) sont des 
contributions pond^r^es ayant pour composantes i (Oii<lst) les 
produits de convolution entre la contribution n ^ la sequence 

15 d' excitation et la r^ponse impulsionnelle b du filtre de 

synthase pond6r6 ; 

- b d^signe le vecteur ligne compost des nc produits 
calaires entre le vecteur X et les vecteurs ligne F^^^^: 

- B d6signe une matrice syin6trique k nc lignes et nc 
colonnes dont le terme B^, j=Fp(i) .Fp( j ^'^ {Oii,j<nc) est 6gal 
au produit scalaire entre les vecteurs Fp^^) et Fp(j) 
pr6c6deinment d6finis ; 

- (.j*^ d6signe la transposition matricielle. 

Pour les impulsions de 1* excitation stochastique dsn 
25 snp=nc-l) les vecteurs Fp^^) sont simplement constitu6s par 
le vecteur de la r6ponse impulsionnelle h d6cal6e de p{n) 
^chantillons. Le fait de tronquer la r^ponse impulsionnelle 
comme d6crit pr6c6demment permet done de r^duire sensiblement 
le nombre d' operations utiles au calcul des produits sca- 
30 laires faisant intervenir ces vecteurs Fp^^,) . Pour la contri- 
bution predictible de 1 ■ excitation, le vecteur Fp(0)=YTP a 
pour composantes Fp^o)^^) (OSi<lst) les convolutions y^p ( i ) 
que le module 38 a calcul6es suivant la formule (1) ou d') 
pour le retard de prediction k long terme select ionn6 TP. Si 
35 MV=0. la contribution n=0 est 6galement de type impulsion- 
nelle et la position p(0) est k calculer. 
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Minimiser 1 ' erreur quadratique E d^finie ci-dessus 
revient trouver 1' ensemble des positions p(n) qui 

maximisent la correlation normalis6e b. B"*^ .b*^ puis ^ calculer 
les gains selon g=b.B"-^. 
5 Mais une recherche exhaustive des positions d' impul- 

sion n6cessiterait un volume de calculs excessif . Pour att6- 
nuer ce problfeme, 1 • approche multi-impulsionnelle applique 
en g6n6ral une procedure sous-optimale consistant h calculer 
successivement les gains et/ou les positions d'irr^ulsion pour 

10 chaque contribution. Pour chaque contribution n {Oin<nc) , on 
determine d'abord la position p(n) qui maximise la correla- 
tion normalis^e (^p^ ^ ^ ^ ^^p' ^p^ ^ ' on recalcule les gains 

^n^0> ^ ^n<^' selon gn^b^.Bj^"^, ou gn^^gn^^^ g^fn)), 

^n~ ^t)(0) , . , . ,b(n) ) et Bj^={Bj^^ j^Oii, jiCn' puis on calcule pour 

15 1' iteration suivante le vecteur-cible e^^ egal au vecteur- 
cible initial X auquel on retranche les contributions 0 A n 
du signal synth6tique pond6re multipli^es par leurs gains 
respectifs : 

n 

A 1' issue de la derni^re iteration nc-1, les gains 
^nc-1^^^ sont les gains seiectionn^s et 1 ' erreur quadratique 
minimisde E est egal ^ I'^nergie du vecteur-cible ©nc-l* 

2 5 La methode ci-dessus donne des r^sultats 

satisf aisants, mais elle n^cessite 1* inversion d'une matrice 
Bj^ A chaque iteration. Dans leur article "Amplitude 
Optimization and Pitch Prediction in Multipulse Coders" (IEEE 
Trans, on Acoustics, Speech, and Signal Processing, Vol.37, 
30 NO 3, Mars 1989, pages 317-327), S. Singhal et B.S, Atal ont 
propose de simplifier le probieme de 1' inversion des matrices 
By^ en utilisant la decomposition de Cholesky : B^^M^^.M^ ou 
Mj^ est une matrice triangulaire inferieure. Cette 
decomposition est possible parce que B^^ est une matrice 

3 5 symetrique A valeurs propres positives. L'avantage de cette 

approche est que 1 ' inversion d'une matrice triangulaire est 
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relativement peu complexe. B„'^ pouvant etre obtenue par 

La decomposition de Cholesky et 1' inversion de la 
matrice Mj^ n^cessitent toutefois d'effectuer des divi-sions 
et des calculs de racines carries qui sont des operations 
exigeantes en termes de complexity de calcul . L' invention 
propose de simplifier consid6rablement la mise en oeuvre de 
1* optimisation en modif iant la decomposition des matrices B 
de la facon suivante : 



n 



ou est une matrice diagonals et L.^ est une matrice 



1 1 

triangulaire inf6rieure n'ayant que des 1 sur sa diagonale 
principale (soit L^^Wj,. K^^/-^ avec les notations 
precedentes) . Compte-tenu de la structure de la matrice Bj^, 
les matrices L^^Rn-^n' ^n' ^n ^n"^ construites 

chacune par simple adjonction d'une ligne aux matrices 
correspondantes de 1' iteration pr^c^dente : 



20 



'I2-1 



B{n,0) \ 



B{n,n-2.) 

^B{n,Q) . . . B{n,n-l) B{n,n) ) 
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n 



[L{n,0) 



-J3-1 



L{,n,n-1) 1 



R(n,n-1) R{n,n)j 
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^12-1 



0 ^ 



^0 . . .0 K{n)j 



10 



.-1 _ 



r -1 



.-1 



Dans ces conditions, la decomposition de B^, 
15 1- inversion de L^^, I'obtention de Bj^"^=iC^. (L^^"^ ) L^^"^ et le 
recalcul des gains ne n^cessitent qu'une seule division par 
iteration et aucun calcul de racine carr6e. 

L* analyse stochastique relative ^ une sous-trame 
d'une trame vois^e (MV=1,2 ou 3) peut d^s lors se d^rouler 
2 0 comme indiQu6 sur les figures 8 ^ 11. Pour calculer le gain 
de prediction k long terme, 1 • index de contribution n est 
initialise ci 0 ^ 1 * etape 180 et le vecteur Fp^Qj est pris 
6gal k la contribution h long terme Y^p fournie par le module 
38, Si n>0, 1 * iteration n commence par la determination 182 
25 de la position p(n) de 1' impulsion n qui maximise la 
quantite : 

min(U2+p, ist) -1 2 
{ 2 hik'p) .e(Jc)) 

30 . P ^ P P ^ 

min{Li2+p, ist) -1 
k-p 

ou e= (e ( 0 e ( lst-1) ) est un vecteur-cible calcuie lors 
35 de 1* iteration precedente. Differentes contraintes peuvent 
§tre apportees au domaine de maximisation de la quantite ci- 



wo 96/21220 



PCT/FR96/00006 



32 

dessus inclus dans I'intervalle [O.lstC L' invent ion utilise 
de pr6f6rence une recherche segnmentaire dans laguelle la 
sous-trame d' excitation est subdivis6e en ns segments de m§me 
longueur (par exemple ns=10 pour 1st =40) . Pour la premiere 
5 impulsion (n=l), la maximisation de (Fp.e*^) 2/ (Fp-Fp"^) est 
effectu6e sur 1' ensemble des positions possibles p dans la 
sous-trame. A 1' iteration n>l, la maximisation est effectuee 
h l'6tape 182 sur 1 ' ensemble des positions possibles k 
1' exclusion des segments dans lesquels ont 6t6 respect ivement 

10 trouv6es les positions pd) p{n-l) des impulsions lors 

des iterations precedentes. 

Dans le cas ou la trame courante a ete detect^e comme 
non vois6e. la contribution n=0 est ^galement constituee par 
une impulsion de position p(0). L'^tape 180 comprend alors 
15 seulement 1 ' initialisation n=0, et elle est suivie par une 
6tape de maximisation identigue a 1 ' 6tape 182 pour trouver 
p(0), avec e=e_i=X comme valeur initiale du vecteur-cible. 

On remarque que lorsque la contribution n=0 est pr^- 
dictible (MV=1. 2 ou 3) . le module 3 8 d' analyse LTP en boucle 
20 ferm^e a effectu6 une operation de nature semblable k la 
maximisation 182, puisqu'il a determine la contribution k 
long terme. caract6ris6e par le retard TP, en maximisant la 
quantity (YT.e'^)2/ (Y,j..Yt'^) dans 1 • interval le de recherche des 
retards T, avec e=e_i=X comme valeur initiale du vecteur- 
25 cible. On peut 6galement, lorsque I'^nergie de la 
contribution LTP est tr^s faible, ignorer cette contribution 
dans le processus de recalcul des gains. 

Apr^s l*6tape 180 ou 182, le module 40 proc^de au 
calcul 184 de la ligne n des matrices L, R et K intervenant 
30 dans la decomposition de la matrice B, ce qui permet de 
completer les matrices L^, Rn et K„ d6finies ci-dessus. La 
decomposition de la matrice B permet d'6crire : 



35 



J-1 

B[n,j) = R{n,j) + X.(2i,Jc) .RU.k) 

Jc=0 



wo 96/21220 



PCT/FR96/00006 



10 



20 



33 

pour la coit^osante situ6e h la ligne n et ^ la colonne j . On 
peut done 6crire, pour j croissant de 0 ^ n-1 : 

Rin.j) 'B{n,j) -J^ L{n.k) .R(j,k) 
Ii(n,J) = R(n,j).K(:J) 

et, pour j=n : 

r ""'^ ^ 

K{n) ^1/R{n,n) =l/[B(n,ja) -J^ L{n,k) .R(n,k)\ 

lj(n,n) « 1 



Ces relations sont exploit^es dans le calcul 184 
d6taill6 sur la figure 9. L* index de colonne j est d'abord 
15 initialise ^ 0, a l*6tape 186. Pour 1 • index de colonne j, la 
variable tmp est d'abord initialis6e a la valeur de la 
composante B(n,j), soit : 



niin(U2+p(j2) .Lh-^pij) ,2st) -1 

X h{k-p{n)) .hik-pij)) 

Jc=max(p(j2) ,p{j) ) 

25 A l'6tape 188, I'entier k est en outre initialise A 

0. On effectue alors une comparaison 190 entre les entiers 
k et j. Si k< j , on ajoute le terme L (n, k) . R ( j , k) ^ la va- 
riable tmp, puis on incr6xnente d'une unit6 I'entier k (etape 
192) avant de r6ex6cuter la comparaison 190, Quand la compa- 

30 raison 190 montre que k= j , on effectue une comparaison 194 
entre les entiers j et n. Si j<n, la composante R(n,j) est 
prise 6gale ci tmp et la composante L(n,j) ^ tmp.K(j) ^ 
1 ' 6tape 196, puis 1 • index de colonne j est increments d'une 
unite avant qu'on revienne ^ l*6tape 188 pour calculer les 

3 5 composantes suivantes. Quand la comparaison 194 montre que 
j=n, la composante K(n) de la ligne n de la matrice K est 
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calcul6e, ce qui termine le calcul 184 relatif k la ligne n. 
K(n> est pris 6gal h l/tmp si tmp#0 (6tape 198) et A 0 sinon. 
on constate que le calcul 184 ne requiert qu'au plus une 
division 198. pour obtenir K(n). En outre, une ^ventuelle 
singularity de la matrice B„ n'entralne pas d' instabilit^s 
puisqu'on 6vite les divisions par 0. 

En r6f6rence k la figure 8. le calcul 184 des lignes 
n de L. R et K est suivi par 1 • inversion 200 de la inatrice 
L constitute des lignes et des colonnes 0 ^ n de la matrice 
l"" Le fait que L soit triangulaire avec des 1 sur sa diago- 
nale principale en simplifie grandement 1- inversion conime le 
montre la figure 10. On peut en effet 6crire : 

L-^'in.j')- -Hn.j')- J: L-^{k',j').L{n.k') 



25 



= -L{n,j')- E L{k',j').L-^{n,k') 
k'^j'*1 



(5) 



pour Q^3'<n et L-l(n.n)=l, c'est-^-dire que 1 • inversion peut 
etre faite sans avoir k op6rer une division. En outre, comne 
les coinposantes de la ligne n de L'^ suffisent k recalculer 
les gains, 1 • utilisation de la relation (5) permet de faxre 
1- inversion sans avoir a m6iaoriser toute la matrice L" , mais 

seulement un vecteur Linv= (Linv(O) Linv{n-1)) avec 

Linv{jM=L-^(n. j") . L' inversion 200 commence alors par une 
initialisation 202 de 1 • index de colonne j' a n-1. A l'6tape 
204, le terme Linv(j') est initialise a-L(n,3M et I'entier 
k' k j'+l. on effectue ensuite une comparaison 206 entre les 
entiers k- et n. Si k-<n, on retranche le terme 
L(k' , j ■ ) .Linv(k' ) ^Linv(j'), puis on incrtmente d ■ une unitt 
I'entier k' (6tape 208) avant de r6ex6cuter la comparaison 
206 Quand la comparaison 206 montre que k'=n, on compare j" 
^ 0 (test 210). Si j'>0, on d6cr6mente I'entier j' d'une 
unit§ (6tape 212) et on revient k 1 • 6tape 204 pour calculer 
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la composante suivante. L' inversion 200 est tennin^e lorsque 
le test 210 montre que j'=0. 

En r6f6rence A la figure 8 1 ' inversion 200 est suivie 
par le calcul 214 des gains r6optimis6s et du vecteur-cible 
5 E pour 1' iteration suivante. Le calcul des gains r6optimis6s 
est 6galeinent trfes simplifi6 par la decomposition retenue 
pour la matrice B. On peut en effet calculer le vecteur 
^n= ^^n^^' ' • • • '^n^^^ * solution de gj^.B^^sb^ selon : 

i=0 

et ' ) =^n-l ' * (n, i ' ) -g^ (n) pour Oii • <n. Le calcul 214 

est d6taill4 sur la figure 11. On calcule d'abord la 
composante b(n) du vecteur b : 

rain(U2+p(ja) ,1st) -1 
^<«> ^^pin) E h{k-p(n) ) .x{k) 

b(n) sert de valeur d' initialisation pour la variable tmq. 

20 A l*6tape 216, on initialise 6galement 1 • index i ^ 0 . On 
effeccue ensuite la comparaison 218 entre les entiers i et 
n. Si i<n, on ajoute le terme b(i) . Linv(i) h la variable tmq 
et on incr6mente i d'une unit6 (6tape 220) avant de revenir 
^ la comparaison 218. Quand la comparaison 218 montre que 

25 i=n, on calcule le gain relatif ^ la contribution n selon 
g (n) =tmq.K(n) , et on initialise la boucle de calcul des 
autres gains et du vecteur-cible (4tape 222) en prenant 
e=X-g(n) .Fpjj^j et i'=:0. Cette boucle comprend une comparaison 
224 entre les entiers i' et n. Si i'<n, le gain g(i') est 

30 recalcul6 ^ 1 ' 6tape 226 en ajoutant Linv (i • ) . g (n) k sa valeur 
calcul^e lors de 1' iteration pr6c6dente n-1, puis on retran- 
che au vecteur-cible e le vecteur 9 (i * ) • Fp ( i . ) . L'6tape 226 
comprend 6galement 1 * incrementation de 1 ' index i* avant de 
revenir ^ la comparaison 224. Le calcul 214 des gains et du 

3 5 vecteur-cible est termini lorsque la comparaison 2 24 montre 
que i*=n. On voit que les gains ont pu Stre mis ^ jour en ne 
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Eaisant appel la ligne n de la ma.rice inverse V 

Le calcul 214 est suivi par une incr&nentacion 228 
de 1 • index n de la contribution, puis par une comparaison 230 
entre 1' index n et le noinbre de contributions nc. S. n<nc. 
5 on revient l> 1 ' 4tape 182 pour !• iteration suivante. 
t optimisation des positions et des .ains est teoarn.e 

lorsaue n=nc au test 230. ^. . 

La recherche segmentaire des itopulsions diminue 
sensiblement le nombre de positions d' impulsion 4 evaluer au 
0 cours des stapes 182 de la recherche de 1 ■ excitation stochas- 
Z... Klle per^et en outre une quantification efficace des 
positions trouvies. Pans le cas typiaue oC la sous-trame de 
lst=40 6=hantillons est divis^e en ns=10 segments de ls=4 
^chantillons. 1 ■ ensemble des positions d' impulsion possibles 
15 peut prendre ns , . Is-P/ Inp. (ns-np, : ) -258 048 valeurs si np=5 
Vl 2 ou 3, =u 860 160 si np=6.MV=0,, au lieu de 
ist!/ npi(lst-np,il=658 008 valeurs si np=5 ou 3 838 380 si 
np=6 dans le cas o^. on impose seulement ,ue deux impulsions 
^ Z puissent pas avoir la m.me position. En d-autres termes 
,0 on peut quantifier les positions sur 18 ^i",- 

bits si np=5, et sur 20 bits au lieu de 22 sr np=6. 

Le cas particulier oi le nombre de segments par sous- 
trame est ^gal au nor^re d' impulsions par 
stochastlque ,ns=np, conduit * la plus grande --P^— ^ ^ 
25 la recherche de 1 ■ excitation stochastique. a nsx ^ - P - 
faible d^bit binaire (si lst=40 et np=5. il y a 8 =3216€ 
ens!illes de positions possibles, quantifiables sur 15 b.t. 
Teulement au lieu de 18 si ns=10, . Mais en r.duisant a c 
point le nombre de sequences d' innovation possibles, on peu. 
30 Lpauvrir la quality du codage. Pour un no-^re d' impulsron, 
ZZ. le no^^re des segments peut *tre optimis* selon u, 
compromis vis^ entre la qualit* du codage et sa simpl.c.f 
de mise en oeuvre lainsi Que le d*bit requis) . 

Le cas oi. ns>np presents en outre I'avantage qu o 
35 peut obtenir une bonne robustesse aux erreurs de transmissio 
en ce qui concerne les positions des impulsions, grace a un 
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quantification s6par4e des num^ros d'ordre des segments 
occup6s et des positions relatives des impulsions dans chaque 
segment occup6 . Pour une impulsion n, le numdro d'ordre s^^ 
du segment et la position relative prj^ sont respect ivement 
5 le quotient et le reste de la division euclidienne de p(n) 
par la longueur Is d'un segment : p (n) =Sj^ . Is+pr^^ (OiSj^<ns, 
Oipry^<ls). Les positions relatives sont chacune quantifi6es 
s^par^ment sur 2 bits, si ls = 4. En cas d*erreur de 
transmission affectant 1 ' un de ces bits, 1 • impulsion 

10 correspondante ne sera que peu d6plac6e, et 1 • impact 
perceptuel de 1 ' erreur sera limits, Les num^ros d ' ordre des 
segments occup^s sont rep6r6s par un mot binaire de ns=10 
bits valant chacun 1 pour les segments occup^s et 0 pour les 
segments dans lesquels 1* excitation stochastique n'a pas 

15 d' impulsion. Les mots binaires possibles sont ceux ayant un 
poids de Hamming de np; ils sont au ncmbre de ns!/tnpi{ns- 
np)!]=252 si np=5, ou 210 si np=6. Ce mot est quantifiable 
par un index de nb bits avec 2'^^~^<ns ! / [np ! (ns-np) ! ) i2^^, 
soit nb=8 dans I'exemple consid6r6. Si, par exemple, 

20 1' analyse stochastique a fourni np=5 impulsions de positions 
4, 12, 21, 34, 38, les positions relatives quantifi^es 
scalairement sont 0,0,1,2,2 et le mot binaire repr^sentant 
les segments occup6s est 0101010011, ou 339 en traduction 
d6cimale. 

2 5 Au niveau du d6codeur, les mots binaires possibles 

sont stock6s dans une table de quantification dans laquelle 
les adresses de lecture sont les index de quantification 
recus . L' ordre dans cette table, d6termin6 une fois pour 
toutes, peut §tre optimist de fagon qu'une erreur de 

3 0 transmission affectant un bit de 1 • index (le cas d* erreur le 

plus frequent, surtout lorsqu'un entrelacement est mis en 
oeuvre dans le codeur canal 22) ait, en moyenne, des 
consequences minimales suivant un crit^re de voisinage. Le 
crit^re de voisinage est par exemple qu'un mot de ns bits ne 
35 puisse Stre remplac6 que par des mots "voisins", 61oign4s 
d'une distance de Hamming au plus egale ^ un seuil np-26, de 
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facon k conserver toutes les impulsions sauf S d'entre elles 
^ des positions valides en cas d'erreur de transmission de 
1. index portant sur un seul bit. D'autres crit^res seraient 
utilisables en substitution ou en complement, par exemple que 
deux mots soient consid6r6s comme voisins si le remplacement 
de 1-un par 1« autre ne modifie pas I'ordre d ■ affectation des 
gains associ6s aux impulsions. 

A des fins d • illustration, on peut considerer le cas 
sin?>lifi6 ou ns=4 et np=2. soit 6 mots binaires possibles 
quant ifiables sur nb=3 bits. Dans ce cas, on peut v6rxfxer 
que la table de quantification pr6sent6e au tableau II permet 
de conserver np-l=l impulsion bien positionnee pour toute 
erreur affectant un bit de 1 • index transmis. II y a 4 cas 
d-erreur (sur un total de 18), pour lesquels on regoit un 
index de quantification qu'on sait &tre erron6 (6 au lieu de 
2 ou 4 • 7 au lieu de 3 ou 5) , mais le d6codeur peut alors 
prendre des mesures limitant la distorsion, par exemple 
r6p6ter la sequence d' innovation relative h la sous-trame 
pr6c6dente ou encore af fecter des mots binaires acceptables 
aux index -impossibles" (par exemple 1001 ou 1010 pour 
1- index 6 et 1100 ou 0110 pour 1 • index 7 conduisent encore 
^ np-l=l impulsion bien positionn6e en cas de r^ceprion de 
6 ou 7 avec une erreur binaire) . 

Dans le cas g6n6ral, 1 ' ordre dans la table de 
quantification des mots peut §tre determine k partir de 
considerations arithm6tiques ou. si cela est insuffisant, en 
simulant sur ordinateur les scenarios d'erreurs (de facon 
exhaustive ou par un 6chantillonnage statistique de type 
Monte-Carlo suivant le nombre de cas d'erreurs possibles). 

Pour s6curiser la transmission de 1 ' index de quanti- 
fication des segments occup6s . on peut en outre tirer parta 
des differentes categories de protection offertes par U 
codeur canal 22, notamment si le critfere de voisinage ne peut 
gtre v6rifie de facon satisf aisante pour tous les ca^ 
d-erreurs possibles affectant un bit de 1' index. Le modul. 
d-ordonnancement 46 peut ainsi mettre dans la cat^gorie d, 
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protection minimale, ou dans la cat6gorie non prot6g6e, un 
certain nombre nx des bits de 1 • index qui, s'ils sont 
affect6s par une erreur de transmission, donnent lieu k un 
mot erron6 mais v6rifiant le critfere de voisinage avec une 
5 probability jug6e satisf aisante, et mettre dans une cat6gorie 
plus prot6g6e les autres bits de 1' index. Cette facon de 
proc^der fait appel ^ un autre ordonnancement des mots dans 
la table de quantification. Get ordonnancement peut 6galement 
§tre optimist au moyen de simulations si on souhaite 
10 maximiser le nombre nx des bits de 1 • index affect^s h la 
cat^gorie la moins protegee. 



index de quantification 


mot d' occupation 


des segments 


decimal 


binaire 


binaire 




decimal 




natural 


naturel 






0 


000 


0011 




3 


1 


GDI 


0101 




5 


2 


010 


1001 




9 


3 


Oil 


1100 




12 


4 


100 


1010 




10 


5 


101 


0110 




6 


(6) 


(110) 


(1001 ou 1010) 




(9 ou 10) 


(7.) 


(111) 


{1100 ou 0110) 




(12 ou 6) 



25 TABLEAU II 

Une possibility est de commencer par constituer une 
liste de mots de ns bits par comptage en code de Gray de 0 
ci 2^^^-!, et d'obtenir la table de quantification ordonn^e en 
3 0 supprimant de cette liste les mots n ' ayant pas un poids de 
Hamming de np. La table ainsi obtenue est telle que deux mots 
cons6cutifs ont une distance de Hamming de np-2. Si les index 
dans cette table ont une representation binaire en code de 
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Gray, toute erreur sur le bit de poids le plus faible fait 
varier 1 ' index de ±1 et entraine done le remplacement du mot 
d- occupation effectif par un mot voisin au sens du seuil np-2 
sur la distance de Hamming, et une erreur sur le i-i^me bit 
5 de poids le plus faible fait aussi varier 1 • index de ±1 avec 
une probability d' environ 2^-^. En placant les nx bits de 
poids faible de 1 • index en code de Gray dans une cat6gorxe 
non prot6g6e, une 6ventuelle erreur de transmission af f ectant 
un de ces bits conduit au remplacement du mot d' occupation 
10 par un mot voisin avec une probability au moins 6gale ^ 
(1+1/2+. . .+l/2"^-^)/nx. Cette probability minimale d^croit 
de 1 a (2/nb) (1-1/2^^) POur nx croissant de 1 k nb . Les 
erreurs aff ectant les nb-nx bits de poids fort de 1 • index 
seront le plus souvent corrig^es grace k la protection que 
15 leur appliQue le codeur canal. La valeur de nx est dans ce 
cas choisie selon un compromis entre la robustesse aux 
erreurs (petites valeurs) et un encombrement r6duit des 
categories protyg^es (grandes valeurs). 

AU niveau du codeur, les mots binaires possibles pour 
20 reprysenter 1- occupation des segments sont rangys en ordre 
croissant dans une table de recherche. Une table d^indexage 
associe a chaque adresse le numyro d' ordre, dans la table de 
quantification stockye au dycodeur, du mot binaire ayant 
cette adresse dans la table de recherche. Dans I'exemple 
25 simplifiy yvoquy ci-dessus, le contenu de la table de 
recherche et de la table d'indexage est donny dans l€ 
tableau III (en valeurs dycimales) . 

La quantification du mot d' occupation des segments 
dyduit des np positions fournies par le module d' analyse 
30 stochastique 40 est effectuye en deux ytapes par le module 
de quantification 44. Une recherche dichotomique est d' abor< 
effectuye dans la table de recherche pour dytermine- 
1. adresse dans cette table du mot k quantifier. L' index d- 
quantification est ensuite obtenu k 1 • adresse dyterminye dan 
3 5 la table d' indexage puis fourni au module 46 d ' ordonnancemen 
des bits. 
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Adresse 


Table de recherche 


Table d'indexage 


0 


3 


0 


1 


5 


1 


2 


6 


5 


3 


9 


2 


4 


10 


4 


5 


12 


3 



TABLEATT TTT 

10 

Le module 44 ef fectue en outre la quantification des 
gains calculus par le module 40. Le gain grpp est par exemple 
quantifie dans I'intervalle tO;l,6], sur 5 bits si MV=1 ou 
2 et sur 6 bits si MV=3 pour tenir compte de la plus grande 

15 importance perceptual le de ce param^tre pour les trames tr^s 
vois^es. Pour le codage des gains associ6s aux impulsions de 
1' excitation stochast ique, on quantifie sur 5 bits la plus 
grande valeur absolue Gs des gains g(l) , . . . ,g(np) , en prenant 
par exemple 32 valeurs de quantification en progression g6o- 

20 m^trique dans I'intervalle [0; 32767], et on quantifie chacun 
des gains relatifs g ( 1 ) /Gs , . . . , g (np) /Gs dans I'intervalle 
[-1;+1], sur 4 bits si MV=1, 2 ou 3, ou sur 5 bits si MV=0. 

Les bits de quantification de Gs sont places dans une 
cat6gorie prot6g6e par le codeur canal 22, de meme que les 

25 bits de poids fort des index de quantification des gains 
relatifs. Les bits de cjuantif ication des gains relatifs sont 
ordonn^s de facon ^ permettre leur affectation aux impulsions 
associ^es appartenant aux segments localises par le mot 
d' occupation. La recherche segmentaire selon 1' invention 

30 permet en outre de prot6ger de mani^re ef ficace les positions 
relatives des impulsions associ^es aux plus grandes valeurs 
de gain. 

Dans le cas ou np=5 et ls=:4, dix bits par sous-trame 
sont n6cessaires pour quantifier les positions relatives des 
35 impulsions dans les segments. On consid^re le cas ou 5 de ces 
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10 bits sont places dans une cat6gorie peu ou pas protegee 
(II) et ou les 5 autres sont places dans une cat^gorie plus 
prot6g6e (IB). La distribution la plus naturelle est de 
placer le bit de poids fort de chaque position relative dans 
la cat^gorie prot6g6e IB. de sorte que les 6ventuelles 
erreurs de transmission affectent plut6t les bits de poids 
fort et ne provoquent done qu'un d6calage d'un echantillon 
pour 1' impulsion correspondante . II est toutefois judicieux, 
pour la quantification des positions relatives, de consid^rer 
les impulsions dans 1 ' ordre d^croissant des valeurs absolues 
des gains associ^s et de placer dans la cat6gorie IB les deux 
bits de quantification de chacune des deux premieres 
positions relatives ainsi que le bit de poids fort de la 
troisi^me. De cette facon, les positions des impulsions sont 
prot6g6es pr6f erentiellement lorsqu'elles sont associ6es h 
des gains importants, ce qui am61iore la quality moyenne 
particuli^rement pour les sous-trames les plus voisees. 

Pour reconstituer les contributions impulsionnelles 
de 1- excitation, le d6codeur 54 localise d'abord les segments 
au moyen du mot d- occupation recu ; il attribue ensuite les 
gains associ6s ; puis il attribue les positions relatives aux 
impulsions sur la base de 1 • ordre d • iit«>ortance des gains. 

on con«>rendra que les diff6rents aspects de 
1- invention d6crits ci-dessus procurent chacun des 
ameliorations propres, et qu'il est done envisageable de les 
mettre en oeuvre ind^pendamment les uns des autres. Leur 
combinaison permet de r6aliser un codeur de performances 
particuli^rement int^ressantes . 

Dans I'exemple de realisation d6crit dans ce qui 
precede, le codeur de parole k 13 kbits/s requiert de 1 • ordre 
de 15 millions d' instructions par seconde (Mips) en virgule 
fixe, on le r^alisera done typiquement en programmant ur. 
processeur de signal num6rique (DSP) du commerce, de meme que 
le decodeur qui ne requiert que de 1 ' ordre de 5 Mips. 
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KEVENDICATIONS 

1. Proc6d6 de codage ^ analyse par synthase d'un 
signal de parole nuin6ris6 en trames successives subdivis^es 
5 en sous-trames comportant: un nombre determine d' 6chanti lions, 
dans lequel on ef f ectue pour chaque trame une analyse par 
prediction lin6aire du signal de parole pour determiner les 
coefficients d'un filtre de synthase k court terme (60) , et 
une analyse en boucle ouverte pour determiner un degre de 

10 voisement de la trame, et on ef fectue pour chaque sous-trame 
au moins une analyse en boucle ferm6e pour determiner une 
sequence d* excitation qui, soumise au filtre de synthase a 
court terme, produit un signal synthetique representatif du 
signal de parole, chaque analyse en boucle feirmee utilisant 

15 la reponse impulsionnelle d'un filtre compose du filtre de 
synthase ^ court terme et d'un filtre de ponderation 
perceptuelle , 

caracterise en ce que, lors de chaque analyse en 
boucle fermee, on utilise ladite reponse impulsionnelle en 

2 0 la tronquant k une longueur de troncature (Lh) au plus egale 
au nombre {1st) d' 6chanti lions par sous-trame et dependant 
de la distribution energetique de ladite reponse et du degre 
de voisement de la trame. 

2. Procede selon la revendication 1, caracterise en 

25 ce que la reponse impulsionnelle du filtre compose est 
calcuiee sur une longueur totale <pst) superieure au nombre 
(1st) d' echanti lions par sous-trame, en ce qu'on determine 
une longueur minimale L.a telle que I'energie de la reponse 
impulsionnelle calcuiee en tronquant ladite reponse a La 

30 echantillons soit au moins egale a une fraction determinee 
(a) de l*energie de la reponse impulsionnelle calcuiee sur 
ladite longueur totale (pst) , et en ce que la longueur de 
troncature (Lh) est egale h la somme de ladite longueur 
minimale La et d'un terme correcteur (A(MV)) dependant du 

35 degre de voisement de la trame si ladite somme est inf erieure 
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au nombre (1st) d' ^chanci lions par sous-trame. 

3. Proc6d6 selon la revendication 2, caract6ris4 en 
ce que ledit terme correcteur (A(MV)) est une fonction 

croissante du degr6 de voisement . 

4. Proc6d6 selon I'une quelconque des revendicat ions 
1^3. caract6ris6 en ce que le filtre de pond^ration 
perceptuelle a une fonction de transfert de la forme 
W(z>=A{z/Yi)/A(z/72) ^/A{z) d^signe la fonction de 

transfert du filtre de synthase a court terme et 7i et Y2 

sont deux coefficients tels que 0<Y2<Yi<l- 

5. Procdd6 selon la revendication 4. caractdrise en 
ce qu-on repr^sente les coefficients du filtre de synthase 
^ court terme par des param^tres de raies spectrales (LSP) . 
en ce qu ■ on quant if ie lesdits param^tres de raies spectrales 
en ce que, pour const ituer le filtre de synthase k court 
terme auquel est soumise la sequence d' excitation relative 
\ au moins une sous-trame d'une trame, on effectue une 
interpolation entre les paramfetres de raies spectrales 
relatifs a ladite trame et ceux relatifs h la trame 
pr6c6dente, et en ce que, pour calculer la r^ponse 
impulsionnelle du filtre compost, on calcule le filtre de 
synthase a court terme sur la base des param^tres de raies 
spectrales quantifies et interpol^s, tandis qu ' on calcule le 
filtre de pond^ration perceptuelle sur la base des paramfetres 
de raies spectrales interpol^s mais non quantifies. 
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